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ABSTRACT 



Context. Magnetohydrodynamic (MHD) instabilities can play an important role in the structure and dynamics of the 
pulsar magnetosphere. 

Aims. We consider the instability caused by differential rotation that is suggested by many theoretical models. 
Methods. Stability is considered by means of a linear analysis within the frame of the force-free MHD. 
Results. We argue that differentially rotating magnetospheres are unstable for any particular geometry of the magnetic 
field and rotation law. The characteristic growth time of instability is of the order of the rotation period. The instability 
can lead to fluctuations of the emission and enhancement of diffusion in the magnetosphere. 
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1. Introduction 

A global structure of the pulsar magnetosphere is the key 
question to answer for understanding the energy outflow to 
the exterior. Likely, the magnetospheres consist of electron- 
positron plasma with some amounts of ions. This plasma 
can affect the radiation produced in the inner region of the 
magnetosphere or at the stellar surface. Therefore, under- 
standing the properties of a magnetosphere is of crucial 
importance for the interpretation of observations. In re- 
cent years, some progress has been achieved in theoretical 
models of the pulsar magnetosphere (see, e.g., Goodwin 
et al. 2004; Contopoulos et al. 1999; Komissarov 2006). 
Apart from a quasi-static structure, however, various non- 
stationary phenomena (such as waves, instabilities, etc.) 
can play an important role. They also may affect the ra- 
diation but, perhaps, idealised quasi-static magnetosphcric 
models cannot be valid in the presence of physical insta- 
bilities. For example, the electrostatic oscillations with a 
low frequency have been considered recently by Mofiz et al. 
(2011), who found that the thermal and magnetic pressures 
can generate oscillations that propagate near the equator. 
These low-frequency electromagnetic waves are of central 
importance for understanding the underlying processes in 
the formation of the radio spectrum (see, e.g., Melrose 1996 
and reference therein). 

Instabilities of magnetohydrodynamic modes can also 
occur in the pulsar magnetosphere. For typical values of the 
magnetic field, the electromagnetic energy density is much 
greater than the kinetic and thermal energy density. This 
suggests that for much of the magnetosphere, the force-free 
condition is a good approximation for determining the mag- 
netic field. Magnetohydrodynamic (MHD) processes under 
this condition are poorly studied but they might be very 
particular. One of the MHD phenomena that can occur in 



the pulsar magnetosphere is the so-called diocotron insta- 
bility, which is the non-neutral plasma analog of the Kelvin- 
Helmholtz instability. This instability has been studied ex- 
tensively in the context of laboratory plasma devices (see, 
e.g., Levy 1965; Davidson 1990; Davidson & Felice 1998). 
The existence around pulsars of a differentially rotating 
equatorial disc with non-vanishing charge density could 
trigger a shearing instability of diocotron modes (Petri et 
al. 2002). In the non-linear regime, the diocotron instabil- 
ity can cause diffusion of the charged particles across the 
magnetic field lines outwards (Petri et al. 2003). The role 
of a diocotron instability in causing drifting subpulses in 
radio pulsar emission has been considered by Fung et al. 
(2006). Note that the diocotron modes should be substan- 
tially suppressed in a neighbourhood of the light cylinder 
where relativistic effects become important (Petri 2007). 

Recently, one more mode of magnetospheric oscillations 
has been considered by Urpin (2011). This mode is closely 
related to the Alfvenenic waves of the standard magnetohy- 
drodynamics modified by the force-free condition and non- 
vanishing electric charge density. Like Alfven waves, the 
magnetospheric waves of a small amplitude are transverse 
(plasma motions are perpendicular to the wavevector). The 
magnetospheric waves can be unstable because there is a 
number of destabilising factors in the magnetosphere (dif- 
ferential rotation, electric currents, non-zero charge density, 
etc.). The electric current usually provides a destabilising 
effect that leads to the so-called Tayler instability (see, e.g., 
Tayler 1973a, b). This instability is well studied in both 
laboratory and stellar conditions. It arises on the Alfven 
time scale and is particularly efficient if the strengths of 
the toroidal and poloidal field components differ substan- 
tially (see, e.g., Bonanno & Urpin 2008a. b). This condition 
can be fulfilled in many magnetospheric models (see, e.g., 
Contopoulos et al. 1999). 
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Various models of the magnetosphere predict that ro- 
tation should be differential (e.g., Mestel & Shibata 1994; 
Contopoulos et al.1999). It is known that differential ro- 
tation in combination with the magnetic field leads to 
the magnetorotational instability (Vclikhov 1959; Balbus 
& Hawley 1991). This instability may occur even in a very 
strong magnetic field (Urpin & Riidiger 2005). However, 
the instability in the pulsar magnetosphere can differ qual- 
itatively from the standard magnetorotational instability 
because of non-vanishing charge density and the force-free 
condition. In this paper, we consider the instability in the 
pulsar magnetosphere caused by differential rotation. 



2. Basic equations 

Despite uncertainties in estimates of many parameters, 
plasma in the pulsar magnetosphere is likely collisional and 
the Coulomb mean free path of particles is shorter than the 
characteristic length scale. Therefore, the MHD description 
can be applied to such plasma. The partial MHD momen- 
tum equations for the electrons and positrons can be ob- 
tained by multiplying the Boltzmann kinetic equation by 
the velocity and integrating over it. Assuming that plasma 
is non-relativistic, the momentum equation for particles of 
the sort a (a — e,p) reads 



the magnetic pressure. Hence, both these forces can be ne- 
glected in Eq.(l) as well. Then, we have 



e a n a E H x B + R 



0. 



(3) 



It is more convenient to use linear combinations of Eq. (3) 
rather than to solve partial equations. Let us define the 
hydrodynamic velocity and electric current as 

V = ^(n e V e + n p Vp), j = e(n p V p - n e V e ), (4) 

where n — n e +n p . Then, the velocities of the electrons and 
positrons are 



1 

2n e 



2n v 



(5) 



If n is much greater than the charge number density 
\n p — n e \, we have V 3> j/en. However, this inequality 
cannot be valid if \n p — n e \ ~ n. In the general case, the 
sum of electron and positron Eqs.(3) yields the equation of 
hydrostatic equilibrium 



p e E + - j x B = 0, 



(6) 



m a n a V a + (V„ • V)V Q = -Vp Q + n a F a + 



e a n a ( E + x B ) + R Q 



(1) 



(see, e.g., Braginskii 1965, where the general plasma formal- 
ism is developed); the dot denotes the partial time deriva- 
tive. Here, V Q is the mean velocity of particles a, n a and p a 
are their number density and pressure, respectively, F a is 
an external force acting on the particles a (in our case F a is 
the gravitational force) , E and B are the electric and mag- 
netic fields, respectively: R a is the internal friction force 
caused by collisions of the particles a with other sorts of 
particles. Since R a is the internal force, the sum of R Q over 
a is zero in accordance with Newton's third law. Hence, we 



R 



. p . 



have in the electron-positron plasma R e 

A calculation of R a is a very complicated problem of 
plasma physics but we will obtain it using simple physical 
arguments. Generally, R Q is proportional to the difference 
of partial velocities of particles (V e — V p ) and to the tem- 
perature gradient (see, e.g., Braginskii 1965). We will ne- 
glect thermal diffusion in R Q because it is usually small in 
astrophysical conditions and take into account only friction 
caused by a difference in the partial velocities. Then, 



Re (^e ^2?) 5 



(2) 



where r e is the relaxation time of electrons. Note that this 
simple expression for R e is often used even in a laboratory 
plasma (Braginskii 1965) and yields qualitatively correct 
results. We assume that accuracy of Eq.(2) is sufficient in 
the pulsar magnetosphere as well. 

In the electron-positron plasma, both sorts of particles 
have a small mass and the inertial force can be neglected 
in Eqs.(l). The gravitation force is also weak because of 
a small mass, and the gas pressure is much smaller than 



where p e = e(n p — n e ) = eSn is the charge density. Taking 
the difference of Eq.(3) for electrons and positrons, we ob- 
tain 



j = Pe V + a[ E + -xB I , 



(7) 



where a = e 2 n p T e /m e is the conductivity of plasma. 
Eqs.(6) and (7) can be written as 



Pe E* + -JxB = 0, J - ctE* = 0, 

c 



where 
E 



V 

EH xB 

c 



J=j-PeV. 



Eliminating E* from Eq.(8) in favour of J, we have 



J H J x B = 0. 

Cpe 



(8) 



(9) 



(10) 



It follows immediately from this equation that J|| = 0. 
Calculating the cross production of Eq.(10) and B, we ob- 
tain J± = as well. Then, second Eq.(8) yields E* = 0. 
Hence, Eqs.(6)-(7) are equivalent to the conditions 



j = p e V , E = x B. 



(11) 



This implies that the force-free condition in combination 
with Ohm's law (Eq.(7)) is equivalent to the condition of 
the frozen-in magnetic field and the presence of only advec- 
tivc currents. Or, in other words, the force-free condition 
and Ohm's law are compatible only if the electric current is 
poorly advective and the magnetic field is frozen-in. Note 
that this statement is valid at any relation between the 
electron and positron number densities. 
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3. Instability in the pulsar magnetosphere 

The set of MHD equations complemented by the Maxwell 
equations reads in the force-free pulsar magnetosphere 



V • E = 47rp e 



VxE = — — , 

c at 



„ „ „ „ „ IdE 4tt. 
VB = 0, VxB=-— + — j, 

c at c 

j^peV, Ew-— xB. 



(12) 



Note one important property of steady state magneto- 
spheres (d/dt — 0). Such magnctospheres can exist only if 
the hydrodynamic velocity is non- vanishing, V^O. Indeed, 
let us assume that V = 0. Then, we have from Eqs.(12) 
that j and E arc equal to zero. If the electric field is zero 
then p e is also vanishing. Since j = the magnetic field 
has a vacuum structure (V • B = 0, V x B = 0), which 
means that the magnetosphere docs not exist at all. The 
conclusion that there should exist hydrodynamic flows in 
the magnetosphere is the intrinsic property of the equa- 
tions of the force-free magnetohydrodynamics and is valid 
at any relation between the electron and positron number 
densities. 

The MHD processes governed by Eq.(12) are very par- 
ticular, which point can be illustrated by considering a 
linear instability. We assune that the electric and mag- 
netic fields are equal to E and B in the unperturbed 
magnetosphere. The unperturbed charge density and ve- 
locity are p e o and V , respectively. For the sake of sim- 
plicity, we assume that motions in the magnetosphere are 
non-relativistic, Vo <C c. Linearising Eqs.(12), we can ob- 
tain the set of equations that describes the behaviour of 
modes with a low amplitude. Small perturbations will be 
indicated by subscript 1. For the sake of simplicity, we 
treat instability of an axisymmetric magnetosphere with 
respect to axisymmetric perturbations. We consider pertur- 
bations with a short wavelength and space-time dependence 
oc exp(zwi— zk-r) where uj and k are the frequency and wave 
vector, respectively; the wave vector k has no (^-component. 
A short wavelength approximation applies if the wavelength 
of perturbations, A = 27r/fc, is short compared to the char- 
acteristic length scale of the magnetosphere, L. Note that, 
generally, the instability criteria for short wavelength per- 
turbations can differ from those for global modes with the 
lengthscale comparable to L. Instability of global modes 
is usually sensitive to details of the global magnetosphcric 
structure and boundary conditions, which are quite uncer- 
tain in the pulsar magnetosphere. In contrast, the insta- 
bility of short wavelength perturbations is entirely deter- 
mined by local characteristics of the magnetosphere, which 
are less uncertain. Note also that the boundary conditions 
and instability of global modes can seriously modify a non- 
linear development of short wavelength perturbations, par- 
ticularly if the global modes grow faster than the short 
wavelength modes. However, in this paper we consider only 
a linear stage of instability. 

As explained above, instability can occur because of ci- 
ther differential rotation or electric currents. The structure 
of a pulsar magnetosphere and its magnetic topology is 
quite uncertain even in the axisymmetric model. Therefore, 
we consider in this paper only instability caused by differ- 
ential rotation and neglect effects associated to electric cur- 



rents. We will show that instability caused by differential 
rotation can arise at any magnetic geometry. 

Substituting the frozen- in condition E = —V x B/c into 
the equation cV x E = —dH/dt (Eq.(12)) and linearising 
the obtained induction equation, we have 



iwBi = Vx(V 1 xB (1 +V x Bi). 



(13) 



If the unperturbed velocity is caused mainly by rotation, 
V = sOe v where s is the cylindrical radius and e v the 
unit vector in the (/^-direction, then 



iwBi - se v (Bi-VO)=iB (k • Vi) -iVi(k • B ). 



(14) 



This is a standard equation that describes perturbations of 
the magnetic field in different types of the magnetorota- 
tional instability (see, e.g., Balbus & Hawley 1991). 

Substituting j from Eq.(ll) into the Maxwell equation 
in the second line of Eq.(12) and linearising it, we obtain 



AnpeoVi 



[i(ck x Bi + wEi) + 4irp el V 



(15) 



Using the linearised frozen-in condition and neglecting 
terms ~ (Vq/c), this expression can be transformed into 



47TC ( O e0 



B x V! 



10 -kxBi-^Vo. (16) 

eO PeO 



4irp, 



The perturbation of the charge density can be calculated 
from the equation p e \ = V ■ Ei/47r. We have with the ac- 
curacy in terms of the lowest order in X/L 



Pel = —[B • (k x Vi) - V • (k x Bi) 

47TC 



(17) 



Substituting Eq.(17) into Eq.(16) and neglecting terms ~ 
V 2 /c 2 , we obtain the second equation, coupling Bi and Vi, 

4ircp e0 Vi +icuB x Vi +iV [B • (kx Vi )] = -ic 2 k x B x . (18) 

Two equations (14) and (18) describe the coupled evolution 
of small perturbations of the velocity and magnetic field. 

We study only MHD modes with uj <C ck because 
one can split electromagnetic and hydromagnetic modes 
in this case. At to <~ ck, a consideration becomes cum- 
bersome since electromagnetic and hydromagnetic modes 
are strongly coupled. Apart from this, MHD effects oper- 
ate basically on a timescale longer than the inverse light 
frequency, (ck) , because V <C c in our model. Therefore, 
one can expect that MHD instability should be particularly 
efficient for modes with a relatively low frequency, oj <C ck. 
That is why we consider such modes first. 

Estimating B\ ~ Vi(kB/u>) from Eq.(14), we obtain 
that the second and third terms on the l.h.s. of Eq.(18) 
are small compared to the term on the r.h.s. by a factor 
uj 2 /c 2 k 2 . Neglecting these terms on the l.h.s., we have 

471-peoVi = -ickxBi. (19) 

Modes turn out to be transverse, k • Vi ~ 0. Subsituting 
Eq.(19) into Eq.(14), we obtain 



iwBi — se v (Bi • VO) 



c(k-Bp) 



k x Bi = 0. 



e0 



(20) 



The dispersion equation can be obtained from Eq.(20) in 
the following way. Calculating a scalar product of Eq.(20) 
and VO, we obtain the expression for (Bi • VO) in terms of 
[Bi • (k x VO)]. Substituting this expression into Eq.(20), 
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we can express after some algebra Bi in terms of [Bi • (k x 
VO)]. Then, a scalar product of the obtained equation and 
(k x VO) yields the dispersion relation in the form 



c 4 fc 2 (k-b) 2 

02 



r 2 (k ■ h) 

S ^[e,.(kxVfl)], (21) 



where m = 47rcp e o/-Bo and b = B /-Bo- 

If rotation is rigid and VO = 0, we obtain the dis- 
persion equation of magnetospheric waves considered by 
Urpin (2011). Note that in our case k • Vo = because 
k has no azimuthal component, whereas Vo corresponds 
to rotation and has only the azimuthal component. If 
|k • b| > fc(fi m |sVO|/c 2 fc 2 ) the dispersion relation for mag- 
netospheric waves reads 



W = ± C (k-b)^-. 



(22) 



Since we assume in our consideration that the frequency of 
magnetohydrodynamic modes should be lower than ck, the 
magnetospheric modes exists if 



O m > c(k-b). 



(23) 



This condition can be satisfied for waves with the wave vec- 
tor almost (but not exactly) perpendicular to the magnetic 
field. If the vector k is almost perpendicular to b, it is con- 
venient to denote the angle between k and b in a meridional 
plane as (ir/2—5). Then, kb = k cos(7r/2 — S) w kSip. Then, 
Eq.(23) is satisfied if S < Q m /ck. 

If rotation is differential and s|VO| > (c 2 fc/O m )(k • b), 
the properties of magnetospheric waves can be quite differ- 
ent. The first term on the r.h.s. of Eq.(21) is always positive 
and cannot lead to instability, but the second term can be 
negative for some k. The instability (u 2 < 0) is possible 
only if the wavevector is almost perpendicular (but not ex- 
actly) to the magnetic field and the scalar product (k • b) is 
small but non- vanishing. Only in this case, the second term 
on the r.h.s. of Eq.(21) can be greater than the first one. 
Let us estimate the range of wave vectors that corresponds 
to unstable perturbations, introducing again the angle be- 
tween k and b as (n/2 — S). Substituting this expression into 
Eq.(21) and estimating [e v ■ (k x VO)] ~ fc|VO|, we obtain 
that the second term on the r.h.s. of Eq.(21) is greater than 
the first one if 



<5< 



s \n m vn\ 

c 2 k 2 



(24) 



The angle <5 turns out to be small, and only perturbations 
with a wave vector almost perpendicular to B can be un- 
stable. 

The instability arises if the second term on the r.h.s. 
of Eq.(21) is positive. Therefore, the necessary condition of 
instability reads 



(k-b) 



[e v ■ (k x VO)] > 0. 



(25) 



Since the sign of TO depends on the charge density, the 
necessary condition is 



(k-b)[e v • (k x VO)] > Oor < 



(26) 



in the region of positive or negative charge density, respec- 
tively. It turns out that the necessary conditions (26) can 
be satisfied by the corresponding choice of the wave vector 



at any VO and b. Indeed, since k is almost perpendicular 
to the magnetic field we can represent it as 



k w ±ke ip x b + 5k, 



(27) 



where 8k is a small component of k parallel (or antiparal- 
lel) to b, k 3> 5k (we neglect terms of the order {Sk/k) 2 ). 
Substituting expression (27) into Eq.(26), we obtain for the 
upper sign with the accuracy in terms linear in 5k 



± k(dk • b)(b • VO) < 0. 



(28) 



Obviously, at any sign of (b • VO), one can choose 8k in 
such a way that condition (28) will be satisfied. Condition 
(26) for the region with a negative charge density can be 
considered by analogy. Hence, the necessary condition of 
instability (25) can be satisfied an any differential rotation 
in the regions of positive and negative charge density. 

The characteristic growth rate can be obtained from 
Eq.(21), using estimate (k • b) ~ kSip. Then, 



— ~ |sVO| 



(29) 



where tq is the growth time of instability caused by differ- 
ential rotation. If differential rotation is sufficiently strong 
and |sVO ~ £1, then the growth time of instability is of 
the order of the rotation period. 



4. Discussion 

We have considered the instability of a pulsar magneto- 
sphere caused by differential rotation. The consideration 
was made using the force-free approximation that should 
be satisfied in the magnetosphere with a high accuracy be- 
cause the electromagnetic energy density is much greater 
than the kinetic and thermal energy density. The main re- 
sult of this study is that the differentially rotating force- 
free magnetosphere is always unstable. This conclusion is 
valid for any particular magnetic topology and rotation law. 
The instability considered in this paper is the representa- 
tive of a wide class of the magnctorotational instabilities 
(see, e.g., Velikhov 1959; Balbus & Hawley 1991; Urpin 
& Riidiger 2005) modified by the force-free condition and 
non-vanishing charge density. The typical growth time of 
the instability can be quite short and can reach the rota- 
tion period in the case of a strong differential rotation with 

|sVO| - n. 

Differential rotation is typical for many models of the 
magnetosphere. For instance, in the axisymmetric model by 
Countopoulos et al. (1999) the angular velocity decreases 
inversely proportional to the cylindrical radius beyond the 
light cylinder and even stronger in front of it. For this rota- 
tion, the growth time of instability should be of the order of 
the rotation period. Numerical simulations by Komissarov 
(2006) showed that within the light cylinder, plasma ro- 
tates differentially basically near the equator and poles. 
Therefore, a strong differential rotation should lead to in- 
stability arising in these regions. However, the situation can 
be quite different near the light cylinder where the insta- 
bility can occur in a much wider region. 

The instability considered can be responsible for fluctu- 
ations of the magnetospheric emission with the character- 
istic timescale ~ 1/lo. Hydrodynamic motions accompany- 
ing the instability can be the reason of turbulent diffusion 
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in the magnetosphere. Note that an influence of the diffu- 
sion coefficients should be strongly anisotropic with a much 
higher enhancement in the direction of the magnetic held 
since the velocity of motions across the held is much slower 
than along it. 

Despite the force-free condition that substantially re- 
duces the number of modes that can exist in the magneto- 
sphere, there are still many destabilising factors that can 
lead to instability. Apart from differential rotation, the elec- 
tric current is likely one more important factor of destabili- 
sation. Note that the topology of the magnetic field can be 
fairly complicated in the magnetosphere, particularly in a 
region close to the neutron star. This may happen because 
the field geometry at the neutron star surface should be 
very complex (see, e.g., Bonanno et al. 2005, 2006; Urpin 
& Gil 2004). Therefore, magnetospheric magnetic config- 
urations can be subject to the so-called Tayler instability 
caused by a distribution of currents. This instability in the 
pulsar magnetosphere will be considered elsewhere. 
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